Entanglement distillation by adiabatic passage in coupled quantum dots 
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Adiabatic passage of two correlated electrons in three coupled quantum dots is shown to provide 
a robust and controlled way of distilling, transporting and detecting spin entanglement, as well as 
of measuring the rate of spin disentanglement. Employing tunable interdot coupling the scheme 
creates, from an unentangled two-electron state, a superposition of spatially separated singlet and 
triplet states. A single measurement of a dot population (charge) collapses the wave function to 
either of these states, realizing entanglement to charge conversion. The scheme is robust, with the 
efficiency close to 100%, for a large range of realistic spectral parameters. 

PACS numbers: 03.67.Mn, 03.67.Hk, 03.67.Lx, 73.63.Kv 
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Creation and detection of spin entanglement is a major 
task for quantum information processing^. A particular 
implementation of the processing relics on electron spins 
in coupled quantum dots, proposed as qubits for quantum 
inverters ]2J and for universal gating in quantum com- 
putation [3j- It has been proposed that entangled two- 
electron spin states in quantum dots can be produced by 
tuned quantum gatesEI, 0, 0| , by filtering through time- 
dependent barriers or by projective measurements 
0)13 • Entanglement is proposed to be detected by cur- 
rent noise measurements |9j. Impressive recent progress 
in coherent control of electronic states in quantum dots 
(nl Ull Hil and spin coherence 0] gives strong im- 
petus to these concepts. 

Here we introduce a scheme for performing spin en- 
tanglement distillation. The remarkable feature of the 
scheme is that, unlike previous proposals, it is also capa- 
ble of entanglement detection, transport, as well as dis- 
entanglement measurement, all in a robust way, without 
the need for fine tuning or precise knowledge of spectral 
or pulse parameters. The scheme is based on our finding 
of a strong correlation between adiabatic passage and 
entanglement: a single adiabatic pulse induces entirely 
different adiabatic passages of different Bell (maximally 
entangled) spin states. We call our scheme, which can 
be realized by current experimental techniques, entan- 
glement distillation by adiabatic passage (EDAP). 

We demonstrate the scheme on two electrons in three 
coupled quantum dots. Starting from an unentangled 
pair, a combination of temporal pulse sequences of in- 
terdot couplings spatially separates singlet and triplet 
states. Entanglement is converted to charge whose de- 
tection uniquely gives a triplet or a singlet. The same 
principle is used for entanglement detection. The scheme 
transports triplets only, leaving singlets, allowing for se- 
lective transport of entangled pairs. Our work is moti- 
vated by quantum control techniques for atomic systems 
such as adiabatic passage and stimulated Raman adia- 
batic passage (STIRAP ) 1151 w hich are increasingly used 
in solid state physics |lfi Il7j. For example, STIRAP 
has been recently proposed to transport single electrons 
[lij , and adiabatic passage to spatially separate electron 
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FIG. 1: Entanglement distillation by adiabatic passage. 
Three quantum dots are coupled via electrode-defined bar- 
riers giving tunnel couplings ti2 and £23- The ground state 
energy of dot 2 is shifted by A. The charge on dots 1 or 3 is 
detected by electrometers. On the right the four figures show 
the scheme at work [the light bulb in d) is an electrometer] 



singlet pairs [19|, in coupled dots. 

We model the physics of two electrons in three coupled 
dots (Fig. ^) by the time-dependent Hubbard Hamilto- 
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with the Fermi creation (ajj and annihilation (a^) op- 
erators for dot i (1, 2, and 3) and spin A =t,4, and 
number operators n.i\ = af x ai\. The confining energies 
Ej do not depend on spin. We take e\ = £3 = 0, while set- 
ting an offset for the middle dot £2 = A. The offset can 
be controlled electrostatically, or it can be fixed within a 
useful spectral range as shown below. We take the on-site 
Coulomb repulsion U^.ii — U to be the same for all dots; 
similarly for the off-site interactions Ui t x-,i+i,X' = V , an d 
zero otherwise. Hopping integrals representing interdot 
couplings are t%j . For our system only t\2 and <23 are 
not zero and depend on time t, so that H = H(t). The 
interdot couplings are modulated by electrostatic gates 
defining interdot barriers. The spectral scales are meVs, 
with t <§C U for realistic systems. In the examples below 
we use generic values of U = 1 mcV, V = or 0.1 meV, 
and maximum hoppings smaller than 0.1 meV. Precise 
values will not be relevant. 

The time dependent spectrum of H, in the presence 
of interdot coupling pulses, is shown in Fig. We 
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take Gaussian pulses of the form tij(t) = t exp(£ 2 /2r 2 ), 
where t$ is the maximum pulse strength and r is the 
dispersion. The overlap between the pulses is taken to 
be 2r, the width of one pulse. There are three weakly 
coupled groups of states. The lowest states with energy 
E « consist of electrons occupying mainly dots 1 and 
3. The highest state, of E U + 2A, is for a double 
occupancy of dot 2. The states relevant for EDAP have 
E rs U, A, and comprise electron singlets and triplets on 
neighboring dots. These states are magnified in Fig. 03. 
To simplify notation we introduce the following labels for 
triplet T and singlet S states on dots i and j (assuming 
i < j), as well as for double occupancy states D: 
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= a+a+\0), \T-i)a = o+o+ |0>, 
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= (l/V2)(a+a+ +a+o+)|0), 


(3) 


\Sh 


= (l/V2)a+a+ -a+a+\0), 
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Here |0) is the vacuum. The triplet states \T$ Z ) are la- 
beled by their spin S z . States \Tq) and \S) are spin en- 
tangled. 

We first summarize EDAP steps and then discuss the 
physics in detail. The scheme is shown in Fig. ^ ( a ) 
Start with two uncoupled electrons occupying neighbor- 
ing dots 1 and 2. (b) Raise slowly the energy of the 
middle dot 2 to A being on the scale of U (this step is 
not necessary if A is built in), (c) Apply an overlapping 
pulse sequence of t\2 and £23 (order not relevant). After 
the pulses fade away, A can be switched back to zero, if 
necessary. The resulting state is with a high probability 
a superposition of a singlet state, spread over dots 1 and 
2, and triplet states, on dots 2 and 3. A detection of (the 
absence of) charge in dot 1, collapses the wave function 
to the singlet (triplet). Mathematically, an initial two- 
electron state *f?(t — 0) localized on dots 1 and 2, is a 
superposition 
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*(0)=a|5)i2 + 6|r )i2 + c|ri)i 2 

After EDAP, the state will be 

¥(00) = o'|S)i2 + 6'|T ) 23 + c / |T 1 ) 23 + d'|T_ 1 ) 23 , (7) 

where the primed coefficients are equal to unprimed up 
to a phase factor. The singlet state returns to the initial 
dots while the triplets are transported to dots 2 and 3. As 
a result entanglement is coupled to charge on dots 1 and 
3. The scheme also works as a noninvasive entanglement 
detector. If the initial state is a singlet, the final state 
is the same (up to a phase). If it is a triplet, the state 
is shifted in space. Charge measurement on dots 1 or 
3, which is a nondcmolition measurement for singlet and 
triplet states in the absence of interdot coupling, sepa- 
rates the two. In general, probabilities of finding, say the 
singlet in a given initial state, |a| 2 can be obtained by 




FIG. 2: (a) Temporal evolution of the two-electron spectrum 
(solid lines) of Hamiltonian H in the presence of two over- 
lapping Gaussian pulses (dashed) of t\2 and £23- The spectra 
are plotted for U = 1 meV, V = 0, and A = 0.8 meV. The 
pulses of £12 (i) and £23 (i) have widths r ^ 500 ps. (b) States 
with S z — relevant for EDAP, from the box in a). There 
is a level repulsion (anticrossing) inside the circles, where the 
passage is rapid. At other two crossings there is no repulsion. 
The horizontal line is the trapped state VPi. (c) Counterintu- 
itive passage scheme for ^1 showing probabilities p of finding 
states I To) 12 and |Tb)23- (d) Passage scheme for | S 1 } 12 showing 
the probabilities p of observing \S)i2 and \Di). 



repeating the measurement on the identically prepared 
state, detecting a degree of entanglement. The scheme 
does not, however, discern the individual triplet states 
\Tq) and \T±i) without an additional single-dot control 
(e.g., spin rotation). Finally, the scheme detects disen- 
tanglement and charge decoherence by observing system- 
atic deviations from the expected final states (e.g., de- 
tecting charge on both 1 and 3). 

To demonstrate the scheme we study the evolution of 
each of the states in the superposition of Eq. © . Con- 
sider triplet states first. It is useful to find the eigenstates 
of H whose energies do not depend on ti2 or i2 3 ; in anal- 
ogy with quantum optics, we call these states trapped. 
There are four two-electron trapped states of H: 

fi = sin<^|T )i2 -cosv3|T ) 23 , (8) 

* 2 = sin<^|Xi) 12 -cosv3|Xi) 23 , (9) 

* 3 - sin(^|T_ 1 ) 12 -cos(^|T_ 1 ) 23 , (10) 

* 4 = [\D)i - \D) 2 + \D) a ] /Vs. (11) 

The mixing angle <p = (p(t) is defined by tan ip = ii 2 /te- 
states ¥1 through ¥ 3 have energy V+ A, while ¥4, which 
is trapped only for A = 0, has energy U. As in STIRAP, 
which is a technique for population transfer via trapped 
states 15]. states ¥1 through ¥3 allow the passage of an 
initial triplet state |T)i2 to |T) 23 , or vice versa. Take 
¥1 a an example. If the initial state is |Tb)i2, it will be 
100% in ¥1 for £23 = 0, when ti2 is slowly turned on 
(ip = tt/2). The state is unaltered until a subsequent 
overlapping pulse of £23 will smoothly move the state to 



= |2b)23) after tx2 vanishes (<p — 0). During the 
passage, no state other than the two triplets is popu- 
lated. The numerical calculation is shown in Fig. 
confirming the qualitative picture. In our context this 
pulse sequence (t±2 before £23) can be called counterintu- 
itive, while the opposite order (£23 before £12) intuitive. 
Transfer through ^1 by counterintuitive sequence is ex- 
tremely robust, independent on spectral parameters, as 
long as adiabatic conditions, to be specified, hold. While 
adiabatic passage via $1 is a nine- level process (there 
are nine S z = basis states), the scheme with ^2 and 
^3, for transporting spin unentangled triplets |Xi) and 
|T_i), is an exact analogue of the three level STIRAP. 
Triplet states can also be transferred through intuitive 
sequencing (not via if A is greater than the inter- 
dot couplings. Such a transfer is less robust, but for our 
scheme it is equally satisfactory as counterintuitive, since 
we need A > £12,^23 to transfer singlet states, as shown 
below. The fourth trapped state, ^4, is a superposition 
of doubly occupied states. Because it cannot be manipu- 
lated with interdot couplings, we call this state globally 
trapped. We will not use this state below. 

Singlet states are not part of the trapped states. If 
the initial state is the singlet 15)12, the above scheme in 
general leads to an arbitrary superposition of eigenstates 
of H for isolated dots. There is, however, a window of 
energy offsets A where the final state will be 115)12, up 
to a phase. Consider states |5*)i2, \D)i, and \D) 2 , with 
average energies A + V, U, and 2 A + U, respectively. If 
we make A on the same scale as U, state \D)2, as well 
as all other eigenstates, will not be easily accessible due 
to spectral separation (Fig. |2K) - We have an effective 
two-level system with Hamiltonian (up to a constant) 

H' = ^(A + V -U)a z +V2t 12 (t)a x , (12) 

where a a are the Pauli matrices. The eigenstates are 

*+ = oos(i?/2)|S)i2 + Bin(<?/2)|i?)i, (13) 
= sin (tf/2)|S)i2-cos (tf/2)|-D)i. (14) 

The mixing angle z? = $(t), restricted to [0, ir], is de- 
fined by tantf = 2y/2ti 2 /(A + V -U). The nature of the 
time evolution of the singlet depends critically on A. In 
resonance, A + V ~ U, the singlet is initially a superpo- 
sition of \&+ and After passage of pulse ti 2 the final 
state will be <I/(oo) = |5)i2 cosa + \D)i sina, where the 
pulse area a = J \/2t\2(t)dt. By fine tuning the pulses 
to a — ir, the final state will be \S)i2- 

The above resonant scheme for singles, though allowing 
fine control, is not robust: it requires both the resonance 
condition and precise knowledge of the pulse area. We in- 
stead explore the large spectral window off the resonance. 
For |A + V — U\ > £12, state |5)i2 will be transported 
back to itself, via (This is analogous to adiabatically 
following of a spin along a magnetic field that rotates 
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FIG. 3: a) Calculated probabilities and electron populations 
after EDAP passage as a function of A, with ^(O) = \S)i2 (a 
and b) and \P (0) = a^a^ |0) (c and d). The thin dotted lines 
in c) are for V = 0.1 meV. Pulses are the same as in Fig. |5] 

along y-axis back and forth adiabatically.) Such a pas- 
sage is very robust. The two-level picture is confirmed by 
the numerical calculation with the full Hamiltonian H in 
Fig. 01. 

Figure shows EDAP results as a function of A, for 
two initial states. For the selected r the initial singlet 
returns to the same state at least 90% of times for A > 
0.6 meV (unless at resonance visible by spikes). This is 
closely mirrored by the charge population of the dots 1 
and 3. Dot 2 has always charge one, except for resonance, 
in which dot 1 can be doubly occupied. The influence of 
off-site Coulomb interaction is seen in Fig. The only 
effect is shifting the resonance from A = U to A = U — V. 

What is the condition on the pulse? Passage of state 
|Z) is adiabatic if \(l\hH(t)\k)\ <C (hcuik) 2 , where k are 
other eigenstates of H(t) and wik are the Bohr frequen- 
cies ji^l • We give rough estimates for the limits on pulse 
dispersion r (switching time), based on the qualitative 
criterion that the smallest relevant Bohr period needs to 
be resolved during the passage. EDAP comprises four 
processes: (i) Adiabatic passage of the triplet state. For 
|A| < t , which can be used for triplet transport, this 
is robust if r ^ h/to- In our scheme A > to and the 
smallest relevant Bohr energy is t^/A. Then r > 17, 
where tl = Ti/^q/t) gives the lower limit, (ii) Adia- 
batic passage of the singlet. This is a two-level scheme 
with states separated by ~ to. Thus r > h/to, which 
is within the range of (ii) and need not be considered 
extra, (iii) Rapid passage of the singlet through the an- 
ticrossing at time t given by t — £i 2 (t) = t23{t) (Fig. 
EJd). The level repulsion is small, since 15*) 12 couples to 
IS 1 ) 23 through spectrally distant states such as |5") 13 and 
I D) 2- Interference in the virtual coupling to these states 
further reduces the anticrossing. One can show that the 
level spacing is V g w 2{P/A)(U - A)/(U + A), vanish- 
ingly small at resonance in the oder 0(P/A). Rapid 
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passage occurs for r' < h/V g , where r' = rV g /t Q is the 
time over which the interdot coupling changes by V g , rel- 
evant for resolving the gap. This gives r < tjj where 
T(,=t l x (A/t )(U+ A)/(U- A) is the upper limit. Fi- 
nally, (iv) EDAP has to be performed within the coher- 
ence time of the system, which is, at low temperatures, 
likely in the nanosecond time scale 0, ^| . Consider- 
ing full coherence, the time limitations on the pulse are 
t~l T T u, which for our model parameters is 100 ps 
to 10 ns. Since the lower limit is given by energy i^/A 
which is on the order of the exchange coupling (J = t\j\J 
in the Hubbard model) for our case of A ss U, the times 
are similar to those used for spin-based quantum com- 
puting 0|. The upper limit tjj increases with decreasing 
\U — A|. The scheme will perform quadratically faster 
for larger couplings. 

To identify numerically the regime of applicability of 
the scheme, we define EDAP efficiency w as 

w = |(*(oo)|S) 12 | 2 + |(*M|T ) 2 3| 2 , (15) 

for a state ^>(t) with the initial condition ^(O) — 
a^aj|j0). This definition is insensitive to the relative 
phase change, and to the relative population of the two 
states. The efficiency is plotted in Fig. 0]as a function 
of A and r for the counterintuitive sequence (intuitive 
shows the same picture except at A close to 0). The 
range of applicability, from 100 ps to 10 ns agrees with 
our analytical estimates for our parameters. The graph 
also shows the predicted increase of applicable r with de- 
creasing \U — A|. It is evident that our scheme is very 
robust, covering large range of spectral values and pulse 
times. The horizontal "cut" at A = U = 1 meV indi- 
cates the resonance oscillations of the |5')i2 — \D)i pair 
for which the efficiency depends on the area of the pulse, 
a, and thus on r. The lower limit on r can be further 
reduced by about a decade (to 50 ps for 98% efficiency 
for our parameters) by decreasing the delay between the 
pulses (not shown here). 

Efficiency w can be measured by performing EDAP 
twice: if the first (distillation) passage results in, say, 
singlet, the second (detection) passage should give ab- 
sence of charge on dot 3, if w « 1. Another interesting 
application of EDAP can be in quantifying the influence 
of a charge probe on the charge itself. Say, use EDAP to 
transport triplets via Wi. Since n 2 1 i'i = at all times 
(in fact, \I/ 1 through ^ 3 are the only eigenstates of H(t) 
that are also eigenstates of 712), a measurement of popu- 
lation on dot 2 should not disturb the state. EDAP effi- 
ciency loss is a measure of the invasiveness of the probe. 

In conclusion, we have proposed a robust and realis- 
tic all-electronic scheme for distilling and detecting two- 
electron spin entanglement in coupled quantum dots. 
The scheme converts entanglement to charge, with close 
to 100% efficiency, by spatially separating singlet and 
triplet states within a single superposition. Because the 
entanglement detection is nondestructive the scheme can 
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FIG. 4: Calculated efficiency w as a function of the pulse 
dispersion r and offset A, for the counterintuitive pulse se- 
quence. The darkest window is for efficiency higher than 98%, 
while the second darkest (red) one is for w > 90%. 



be used repetitively to transport entangled pairs as well 
as to detect disentanglement. Although we have used 
the example of quantum dots, wc believe that EDAP is 
general enough to be applicable in other physical imple- 
mentations of quantum information processing. 

This work was supported by the US ONR and FWF. 
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